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ABSTRACT 

This paper demonstrates the feasibility of a procedure for the remote measurement 
of sea-surface temperature which inherently corrects for the effect of the intervening 
atmosphere without recourse to climatological data. The procedure relies upon the near- 
linear differential absorption properties of the infrared window region between 10 and 
13 nm and requires radiometric measurements in a minimum of two spectral intervals 
within the infrared window which have a significant difference in absorption coefficient. 
The procedure has been applied to Nimbus 4 infrared interferometer spectrometer (IRIS) 
data and to Skylab EREP S191 spectrometer data, and it is demonstrated that atmos- 
pheric effects on the observed brightness temperature can be reduced to less than 
1.0 Kelvin. 


INTRODUCTION 

Global measurement of sea-surface temperature (hereafter referred to as SST) on a 
daily basis is currently operational at NOAA's National Environmental Satellite Service 
(NESS) Center. SST data are archived after a rather complex analysis of data from the 
scanning radiometer (SR) aboard the NOAA-4 satellite. The procedure for extracting 
SST information from SR data is an extension of that of Smith and Rao [1]. Basically, 
infrared window measurements (10.5 - 12.5 ^m) are converted to SST values, at de- 
graded spatial resolution, with proper consideration of SR instrument noise, noise re- 
lated to signal transmission, receiving, and processing, and after accounting for the 
effects of the intervening atmosphere based upon climatological data. 

Because of inherent limitations in the single- channel scanning radiometer method 
of producing SST’s, archived temperatures have limited accuracy ^d spatial resolution. 
The next generation of TIROS satellites, denoted TIROS-N series, will contain a five- 
channel advanced very high resolution radiometer (AVHRR), which will provide for im- 
proved' SST retrieval. The first instrument to be flown, planned for FY 78, will contain 
four channels, with a fifth channel to be added on the fourth or fifth' satellite in the 
TIROS-N series. The nominal spectral regions for each of the five channels are 0.55 to 
0.9 /im, 0.72 /im to detector cutoff (approximately 1.0 /im), 3.55 to 3.93 jum, 10.5 to 
11.5 /im, and 11.5 to 12.5 /xm, which will be the fifth channel to be added. The spectral 
response of the first four channels has been finalized, however, the spectral response 
of the fifth channel is still open. Of major impact to SST measurement is the addition 
of the infrared channel from 3.55 - 3.93 nm and the division of the present infrared 
window channel, i. e. , 10.5 - 12.5 /xm, into two channels. The reason for the addition 
of the two infrared channels is to provide additional information to correct for the effects 
of atmospheric moisture on the observed radiance (which can be significant, particularly 
for warm and moist conditions) without recourse to climatological data. 

Previous studies [2,3,4, 5] have indicated that the addition of a second channel in. the 
infrared window, one somewhat less transparent than the other, would allow for compen- 
sation of the effects of absorption and emission by atmospheric water vapor without re- 
course to climatological or other supportive data. The investigations involved the use 
of radiative transfer models, and large field-of-view space acquired data, i. e. , the 
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infrared interferometer spectrometer (IRIS) aboard Nimbus 4. Although the results 
were encouraging, they remain tentative until they can be verified by an experiment for 
which the sensor field-of-view is comparable to that which will be used on future TIROS 
satellites. The Skylab mission provided the first opportunity for such a verification 
experiment. 

This paper describes the results of an investigation which utilized IRIS and SI 91 
spectrometer data acquired during two of the three Skylab missions to validate the 
radiative transfer models used in the initial investigations, and test the validity and 
correctness of two-channel temperature estimating algorithms developed from the radi- 
ative transfer models. The results of the investigation indicated that a significant 
benefit will be derived by the addition of a second channel in the infrared window region. 
Specifically, without recourse to climatological data, the investigation indicated SST's 
could be estimated to within 1°K.^ Although the spectral bands selected for the analysis 
yielded a satisfactory result, the investigation also indicated that the infrared window 
transmission function contains nonlinearities and uncertain aerosol effects which could 
have a significant impact upon the optimum choice of the two spectral channels. Since 
the present analysis was limited in scope it is recommended that a further analysis be 
performed before a final choice of spectral channels is made. 

THEORETICAL BASIS OF REMOTE SST 
MEASUREMENT AT THERMAL INFRARED WAVELENGTHS 

The spectral radiance emitted by an opaque body at wavelength X is given by 

L(X,T) = e(X)L^V,T) (1) 

where c(X) is the spectral emittance of the opaque body and l“°(X, T) is the spectral 
radiance emitted by a blackbody. The latter is represented as 

(3) 

where 

T = the temperature of a blackbody 
c = the velocity of light 
h = Planck's constant 
X = wavelength 
k = Boltzmann's constant 

It is clear from these expressions that if the emittance is known, the temperature 
can be determined by measuring the emitted spectral radiance and inverting Eq. (1). 

The application of such a measurement procedure to determine the temperature of a 
water surface exposed to the atmosphere is more complex. The spectral emittance of a 
sea surface is less than 1 depending upon the emittance angle. Consequently, as one 
attempts to measure the emitted radiation, some sky radiation wiU be reflected from 
tha water rorface and collected by the infrared sensor. Also, since water does not be- 
come opaque to infrared radiation at thermal wavelengths until a depth of approximately 
0.10 mm, some of the measured radiation emanates from below the surface, which 
generally has a slightly different temperature. Therefore, the temperature derived 
from a measurement of the radiance at the surface will be the temperature of a blackbody 
which yields an equivalent value of radiance (i. e. , the "equivalent radiometric tempera- 
ture"). H will be different from the actual surface ten^)erature; the degree of difference 

Temperature accuracies are based upon atmospheric effects only. Other sources 
of noise and accuracy degradation factors are not included. 
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will depend upon the magnitude of the reflected radiation and the temperature gradient 
near the surface. 

The present analysis is not concerned with the relationship between the equivalent 
radiometric temperature and the actual surface temperature, but only with the effect 
of the atmosphere on the equivalent radiometric temperature derived from a radiometric 
measurement performed at satellite altitudes. Therefore, all future references to SST 
will refer to the equivalent radiometric temperature that would be derived from a radi- 
ance measurement at the surface. 

Before reaching a spaceborne sensor, the spectral radiance emanating from the sea 
surface will be atten\iat^ by atmospheric constituents, such as clouds, haze, and ab- 
sorbing gases. These atmospheric constituents also emit and scatter radiant energy, 
which contributes to the total signal received by the sensor. The central problem in 
accurately measuring SST from space lies in determining the extent to which such 
effects can be observed and compensated. 

To demonstrate more clearly the nature of the problem of measuring the sea temper- 
ature from space, consider the spectral radiance leaving the top of the atmosphere, 

L(X), which can be represented by: 

L(X) = T(P )]t(X, P )€(X) + J* L^^[x, T(P)]dr(X, P) (3) 

where 

Po is surface pressure 
T is atmospheric transmittance 
P is atmospheric pressure 
T is temperature 

According to Eq. (3) estimating T(P ) from L(X) requires values of temperature, pres- 
sure, and the differential absorption properties of the atmosphere. Estimating SST 
from a single channel radiometric measurement is analogous to inverting Eq. (3) given 
temperature and pressure values consistent with the prevailing conditions which are ob- 
tained from either atmospheric soundings or climatological data. 

To eliminate the requirement for ancillary data a scheme was devised which utilizes 
more than one spectral channel. The scheme was originally discussed and presented by 
Anding and Kauth [2], and subsequently discussed by McMillin [3] and Prabhakara [4]. 

The technique, like most remote sounding methods, is based on the use of the differential 
optical properties of the atmosphere in the infrared window region to infer the atmos- 
pheric attenuation. The attenuation values are then used to correct for the effect of the 
atmosphere on radiometric data. 

Theoretical Basis of the Multi-Channel Method 

Following the development of Prabhakara [4] the radiative transfer equation may be 
simplified as: 

L(X) = T(P„)]t(X, P^)e(X) + 1.%) [1 - r(X, P„)] (4) 

where L (X) is the weighted mean Planck emission of the atmosphere. 

In the infrared window region there are three primary contributors to the absorption. 
Local water vapor lines, H 2 O continuum, and aerosols. The transmission functions for 
continuum absorption and aerosols are accurately represented by Beers Law. This is 
also true for selective line absorption when the absorption is either weak or the 
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individual lines are heavily overlapped. These conditions are approximately satisfied in 
the window region for one airmass. Hence, window transmission can be expressed as: 

-(k +k +k )u -k.u 
^ ^ -L c a „ t 

T = e = e (5) 

where kt, k^, and k^ are the local line, H 2 O continuum, and aerosol extinction co- 
efficients, respectively, u is the effective absorber thickness. To a good approxima- 
tion Eq. (5) can be represented by the first two terms of its series expansion, i.- e. , 

T 1 - k(X)u (6) 

Substituting (6) into (4), and letting €(X) = 1, we have 


L(X) - L^V»T(Pq)] - |l‘’V.T(Pq)] - L^^(X) I k(x)u . 


(7) 


Expanding the Planck function about the surface temperature T(Pq) and retaining only 
the linear term we have 

bb bb bL^^[X,T(P )] 

T] = T(P^)] + [T - T(Pq)] . (8) 


t)T 

This approximation holds for a small wavelength region and a small range of tempera- 
tures. This relationship allows Eq. (7) to be expressed as: 


T(X) = T(P^) - [T(P^) - T(X)]k(x)u 


(9) 


where T(X) is the brightness temperature of the observed radiance, T(X) is the equiva- 
lent brightness temperature of the atmosphere, and T(Pq) is the surface temperature. 
Eq. (9) shows a lin^r relationship between brightness temperature and absorption co- 
efficient, provided T(X) is not strongly dependent upon X over the spectral region. A 
study by McMillin [3], and results of the present study, show variations of less than 5 
per cent. Conceptually, therefore, measurements in only two wavelength intervals for 
which the respective absorption coefficients are significantly different are required to 
define the linear relationship. ' 1 

To demonstrate the concept three wavelength intervals within the infrared window 
between 10 and 13 were chosen. Three intervals were chosen, rather than only 
two, to observe possible nonlinearities between brightness temperature and absorption 
coefficient. The intervals chosen were those used by Prabhakara [4], and are re- 
spectively, 10.25 - 11.25, 11.25 - 12.0, and 12.0 - 1^9 pm. This result is a conse- 
quence of selecting a 1 pm wide interval beginning at the long wavelength side of the 
9.6 pm ozone band (the most transparent part of the window region) and then dividing 
the remaining wavelength interval into two equal wavenumber segments, each 56 cm“^ 
wide. 

Absorption coefficient evaluation . - The mean absorption coefficient for each 
interval was evaluated empirically from transmission spectra computed from an an- 
aljrtical transmission code. For the infrared window region the code utilizes two 
absorption models; one for H 2 O local line absorption and one for H„0 continuum 
absorption. 

Local line absorption is represented by the Goody model [6] which is given by: 


t(x) = exp 


y; 


-(S/d)u 


2 / S X 


( 10 ) 
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where S/d 
u 
P 
S . 

2770 ! 


intensity to line-spacing parameter (cm”^) 
absorber thickness (pr. cm. ) 

Curtis - Godson equivalent pressure (atm) 

Intensity to half -width parameter (atni . cm”^) 


The parameters S/d and S/27ra were evaluated from a tabulation of spectral line param- 
eters [7] using a procedure discussed by Goody [6], modified to account for an instrument 
slit function. The parameters were ev^uated at a spectral resolution of 10 cm"^ de- 
fined by the width of the slit function when the transmission is 50%. A comparison be- 
tween band model and line-by-line derived spectra is shown in Figure 1. Observe that 
the transmission spectrum is well represented by the Goody model. 


Continuum absorption in the window region results from two mechanisms; that 
caused by the wings of water vapor lines within the 6.3 pm band and the rotational water 
band which are pressure broadened by foreign gases, and that caused by the same water 
vapor lines which are self-broadened. The continuum absorption coefficient at total 
pressure P, and water vapor partial pressure p, is given by: 

k(P,p) =kjP + k2P (11) 


where 

kl is the absorption coefficient for foreign broadening at unit total pressure; 

k 2 is the coefficient for self -broadening at imit water vapor partial pressure. 

The values of kj • and k 2 adapted for the continuum model are based upon a sub- 
jective analysis of the data of Bignell [8], Burch [9], and McCoy [10]. The values 
adopted for k 2 are illustrated in Figure 2. Discussions with Burch and Long indicated 
that a least-squares fit to the data would yield a self -broadening coefficient that was too 
high because of systematic errors for the larger data points and because of the consis- 
tency of the results of McCoy for the CO 2 lasing line at 10.59 pm. Therefore the 
selected values coincide with the lowest values shown. 

The values adopted for k^ are based upon the 10.59 CO 2 laser measurements of 
McCoy [10]. McCoy measured the transmission as a function of total pressure for a 
fixed water vapor pressure. Using the self -broadening coefficient (k 2 ) as a basis, the 
foreign-broadening coefficient (kj) was determined to be 0.005 k 2 . 

Both kj and k 2 are temperature dependent, the value of kj increasing with in- 
creasing temperature and the value of k 2 decreasing with increasing temperature. For 
self -induced absorption the dependence was adopted from the work of Bignell [8] at 2 
percent per degree Kelvin. For foreign-induced absorption the temperature dependence 
was scaled from the temperature dependence of the rotational water lines at 2 percent 
per degree Kelvin. 

The empirical evaluation of the absorption coefficients was accomplished as follows. 
Ten atmospheric representations of temperature and humidity were selected from the 
NASA Four-Dimensional Model Atmosphere compilation [11], representing moist and dry 
conditions (+ one standard deviation from mean moisture conditions) for five oceanic 
global regions extending from the, North Sea to the equator. These atmospheric data were 
input to the transmission code and spectral transmission from 10 to 13 ^m was com- 
puted. These values were integrated over the respective wavelength intervals and trans- 
mission versus effective absorber thickness was plotted. The result is presented in 
F^ure 3. Assuming transmission is accurately represented by Beers Law, i. e. , 

’ T = exp (-ku), linear least squares fits to -tnr versus u yielded three respective 
values for the absorption coefficient, which are also given in Figure 3. 
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To demonstrate the relationship between brightness temperature and absorption co- 
efficient, the ten model atmospheres were input to atmospheric radiance codes and the 
vertical upwelling radiance at space altitudes, emanating from a sea surface whose 
temperature equaled T(P^) was calculated. To establish realism to the computations, 
and to add an additive source of noise, aerosol effects were included in the calculations. 
Specifically, two aerosol representations were selected, taken from the work of Fenn[12]. 
One whose optical properties and size distribution were consistent with a 100 percent 
maritime haze and the other with a 60 percent maritime haze and a 40 percent continental 
haze. The altitude distribution of aerosol density was assumed variable and controlled by 
the sea-level visibility. A 23 km visibility was used for the maritime haze, a 10 km 
visibility for the maritime-continental haze. The calculations of radiance were made by 
numerically evaluating Eq. (3), using the infrared window transmission described above 
to evaluate drCr, P). For each of the twenty spectra, in-band brightness temperatures 
were evaluated and plotted versus absorption coefficient. The results are shown in 
Figure 4a - 4e. The straight lines are least squares fits to the data. 

Discussion of results. - Observe that a nonlinearity exists between brightness tem- 
perature and absorption coefficient, which diminishes at higher brightness temperatures. 
This occurs because the atmospheric brightness temperature increases approximately 5 
percent at the longer wavelength regions, and the increase is observable at the lower 
brightness temperatures, but becomes nearly insignificant for brightness temperatures 
above 290° K. Also note that the atmosphere causes a decrease in the observed bright- 
ness temperature of the sea surface ranging from 2° K for the cool northern region to 
greater than 5° K for the warm equatorial region. The application of the technique 
does, however, estimate the input SST values to within ± 0.5° K. Furthermore, the 
3 -band estimate is insignificantly better than the 2 -band estimate (± 0.5° K compared to 
± 0.3° K). 

The application of the technique to actual space-acquired data is shown in Figures 5 
and 6. For each of the six cases the brightness temperatures were taken from Prabhakara 
[4], which were measured by the IRIS instrument on Nimbus 4. Observe nearly perfect 
linearity is demonstrated and that excellent agreement is obtained between predichons and 
ship measurements for three of the cases (Figure 5). The reason for the discrepancy for 
the other three cases (Figure 6) is unknown, although it could be attributable to a de- 
crease in surface emissivity resulting from a high sea state. 

The utilization of Skylab EREP S191 data as a verification e^eriment fell far short 
of e}q)ectations because of instrument problems. Although the problems have been identi- 
fied, their effects remain uncorrected. Of the eight test areas planned for usage only one 
gave results which were considered satisfactory. This was for the Monroe Reservoir on 
10 June 1973. A circumstance of the spectrometer was that whenever the brightness 
temperature of the spectral radiance was significantly different than either the ambient 
calibration source temperature, or the temperature of the internal mirrors or dichroic, 
inaccurate data resulted. For the Monroe Reservoir the aperture brightness temperature 
was approximately equal to the internal instrument temperatures and accurate data were 
obtain^. 

During the overpass of EREP the surface temperature was measured with a Barnes 
PRT5 radiometer yielding an average value of 25° C. Also, temperature and humidity 
data were obtained by a local radiosonde. Both the surface temperature and radiosonde 
data were input to the radiative transfer model and the radiance at the SI 91 aperture 
was cal ciliated. Also included in the calculation was a 23 km visibility continental haze. 

A comparison of measured and calculated radiance is presented in Figure 7. Generally, 
the agreement is very good. Minor differences are noted between 8 and 9 /im, probably 
caused by the reduced values of aperture brightness temperature. Of major concern for 
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the present study is the region between approximately 10.5 and 13.0 ^m, where the agree- 
ment is within a few percent. It would appear from this comparison that the radiative 
transfer model is a reasonable representation of reality. 

To examine the application of the SST estimation technique to this data, for both the 
EREP spectra and the model calculations, the inband brightness temperatures were com- 
puted and plotted versus absorption coefficient. Least-squares lines were fitted to the 
data and the results are shown in Figure 8. Observe that both the model calculations and 
the EREP data estimate a surface temperature within ± 1.0° K of the recorded ground 
truth temperature. 

CONCLUSIONS AND RECOMMENDATIONS 

The results presented herein indicate that a significant benefit will be derived by the 
addition of a second channel in the infrared window region. Specifically, atmospheric - 
effect uncertainties can probably be reduced to less than 1.0° K without recourse to 
climatological data. The technique is likely to produce a significant benefit for unusual 
conditions, such as warm moist atmospheres over cool waters or vice versa, when 
climatological data would yield a particularly poor result. 

The results presented herein are based upon a study of limited scope, with recourse 
to a limited amount of data. Although the study results are indicative, it would be desir- 
able to apply the technique to a broader data base before a final choice is made for the 
spectral bands to be used in the latter satellites of the TIROS-N series. 

f 
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Figure 3. Average Transmission Versus Effective Absorber Thickness 
for Three Wavelength Regions (10.25 - 11.25, 11.25 - 12.0, 
12.0 - 12.9 /im). ki, k 2 , k 3 are the respective Beers Law 
Absorption Coefficients. 
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Figure 4a. Brightness Temperature Versus Absorption Coefficient. 
Key: Q Low Humidity, Continental-Maritime Hazep; 

® Low Humidity, Maritime Haze; 

A High Humidity, Continental- Maritime Haze;^ 

• High Humidity, Maritime Haze. ‘ 
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Figure 4b. Brightness Temperature Versus Absorption Coefficient. 
Key: ® Low Humidity, Continental-Maritime Haze; 

® Low Humidity, Maritime Haze; 

A High Humidity, Continental-Maritime Haze; 

• High Humidity, Maritime Haze. 
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Figure 4c. Brightness Temperature Versus Absorption Coefficient. 
Key: ® Low Humidity, Continental-Maritime Haze; 

® Low Humidity, Maritime Haze; 

A High Humidity, Continental-Maritime Haze; 

• High Humidity, Maritime Haze. 
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RELATIVE ABSORPTION COEFFICIENT 

Figure 4d. Brightness Temperature Versus Absorption Coefficient. 
Key: ® Low Humidity, Continental-Maritime Haze; 

® Low Humidity, Maritime Haze; 

A High Humidity, Continental-Maritime Haze; 

. High Humidity, Maritime Haze. 
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Figure 4e. Brightness Temperature Versus Absorption Coefficient, 
Key: a Low Humidity, Continental-Maritime Haze; 

0 ;Low Humidity, Maritime Haze; 

4 High Humidity, Continental-Maritime Haze; 
High Humidity, Maritime Haze. 
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Figure 5. IRIS Measured Brightness Temperature 
Versus Relative Absorption Coefficient. 
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Figure 6. IRIS Measured Brightness Temperature 
Versus Absorption Coefficient. 
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Figure 7. Comparison of Measured and Predicted Spectral Radiance for Monroe Reservoir, 
Sa^m, Illinois, on 10 June 1973.. Water Surface Temperature = 298° K. 




Figure 8. Brightness Temperature Versus Absorption Coefficient 
for Monroe Reservoir on 10 June 1973. Model Pre- 
dictions Compared to S191 Measurement Predictions. 
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